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Abstract 

We calculate the leading logarithmic QCD corrections to the matrix element of the 
decay b — > de^e~ in the two Higgs doublet model with tree level flavor changing currents 
(model III). We continue studying the differential branching ratio and the CP violating 
asymmetry for the exclusive decays B — > TTe~^e~ and B — > pe~^e~ and analysing the de- 
pendencies of these quantities on the selected model III parameters, including the 
leading logarithmic QCD corrections. Further, we present the forward-backward asym- 
metry of dileptons for the decay B — > pe~^e~ and discuss the dependencies to the model 
III parameters. We observe that there is a possibility to enhance the branching ratios 
and suppress the CP violating effects for both decays in the framework of the model III. 
Therefore, the measurements of these quantities will be an efficient tool to search the new 
physics beyond the SM. 
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1 Introduction 



Rare B meson decays, induced by flavor changing neutral current (FCNC) b s{d) transitions 
are one of the interesting research area to test the Standard model (SM) at loop level. They are 
informative in the determination of the fundamental parameters, such as Cabbibo-Kobayashi- 
Maskawa (CKM) matrix elements, leptonic decay constants, etc. and useful for establishing the 
physics beyond the SM, such as two Higgs Doublet model (2HDM), Minimal Supersymmetric 
extension of the SM (MSSM) 0, etc. 

Since the SM predicts the large Branching ratio [Br), which is measurable in the near 
future, the exclusive decays induced by 6 — > sl^l^ process become attractive. Such transitions 
has been investigated extensively in the SM, 2HDM and MSSM, in the literature 0- |T^. For 



these transitions, the matrix element contains a term includes the virtual effects of the top 
quark proportional to VtbV^* and additional terms describing the cc and uu loops, proportional 
to VcbV*^ and VubV*^ respectively. Using the unitarity of CKM matrix, i.e. VibV*. = 0, i = u,c, t, 
and neglecting the factor VubV*^ compared to VtbV^* and VcbV*^, it is easy to see that the matrix 
element involves only one independent CKM factor, Vf^V^*. This causes that the CP violating 
effects are suppressed within the SM ||T6| , p!7| . However, for b dl^l~ decay, all the CKM 
factors VtbV^% VcbV*^ and VubV*^ are at the same order and this leads to a considerable CP 
violating asymmetry between the channels induced by the inclusive b —>■ dl^l~ and b dl^l~ 
decays. These effects have been studied in the literature for the inclusive b — > de^e~ decay, in 
the framework of the SM |jl8|. The difficulties of the experimental investigation of the inclusive 
decays stimulate the study of the exclusive decays. However, the theoretical analysis of the 
exclusive decays is complicated due to the hadronic form factors which can be calculated using 
non-perturbative methods. The dispersion formulation of the light-cone constituent quark 
model is one of the method which can be used to calculate the hadronic matrix elements. In 
the literature, the form factors for b de^e~ induced exclusive B (7r,p)e"^e~ decays have 



been calculated in the framework of this method IT9, 201. The CP violation effects for these 



exclusive decays have been studied in the framework of the SM |21 



In this work, we present the leading logarithmic (LLog) QCD corrected effective Hamiltonian 
in the 2HDM with flavor changing neutral currents (model HI) for the inclusive b — > de^e~ 
decay and calculate the differential Br of the exclusive B —>■ (7r,p)e+e^ process. Further, we 
study the CP-violation asymmetry (Acp) and forward-backward asymmetry {A^b) of dileptons 
for the decay B pe~^e~ . 

The paper is organized as follows: In Section 2, we give the LLog QCD corrected Hamilto- 
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nian responsible for the inclusive b de^e~ decay and calculate the matrix element. In section 
3, we present the Br and Acp of the exclusive B — > 7re+e~ decay and analyse the dependencies 
of the Br and Aqp on the couplings In section 4, we study the Br ,Acp and ApB of 

the exclusive B pe^e^ decay . Section 5 is devoted to our conclusions. In Appendix, we 
summarize the essential points of the model III and give the explicit forms of some functions 
we use in our calculations. 

2 Leading logarithmic improved short-distance contri- 
butions in the model III for the decay b de^e~ with 
additional long-distance effects 

In this section, we present the LLog QCD corrections to the inclusive b —>■ de^e~ decay ampli- 
tude in the 2HDM with tree level neutral currents (model III). The LLog QCD corrections to 
the b — i> de^e~ decay amplitude can be calculated using the effective theory. In this method, 
the heavy degrees of freedom, t quark, W^,H^,Hi, and H2 bosons, in the present case, are 
integrated out. Here if^ denote charged. Hi and H2 denote neutral Higgs bosons. The proce- 
dure is to match the full theory with the effective low energy theory at the high scale /i = miy 
and evaluate the Wilson coefficients from mw down to the lower scale /i ~ 0{mi,). In our 
calculations we choose the higher scale as p = my/ since the current theoretical restrictions 
T\ E3] show that the charged Higgs mass is enough heavy to neglect the running from 



to niw 

The effective Hamiltonian relevant for the decay b — > de^e' in the model HI is 

K// = -4%^^,^;, { E (a(/i)o.(/x) + c:(/i)o:(/x)) 

V2 i=l,...,12 

+ \u y: (c^.(/^)(o.(/i)-o-(/x))+c:(^)(o:(^)-o-(^))} (1) 

2=1,2,11,12 

where 0\ ^ , O"'' ■* , are the operators given in eqs. (D, (I) and Cp are the Wilson coefficients 
renormalized at the scale p. Here the unitarity of the Cobayashi-Maskawa matrix (CKM) is 
used, i.e. VtbV^^ + VubV*^ = —VcbV*^ and the parameter A„ is defined as: 

VtbV;, 

Using Wolfenstein parametrization [^, can be written as 



_ p(l-p)-r/^ 7] , 

" (l-p)2+,^2 ^(l_p)2+^2 
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where p, rj and A ~ 0.221 are Wolfenstein parameters. The parameter t] (and therefore A^) is 
the reason for the CP violation in the SM. 

The operator basis is similar to the one used for model III (p5[ and references therein), which 
is obtained by replacing s-quark by rf-quark and adding new operators, i.e. O"*-'-* , z = 1, 2, 11, 12: 

O2 = {dLalt,CLa){cLpl^hLp), 

02 = {dLalfiULa){uLf3Y^Lp), 

03 = idLalfMbLa) i^LfSj^ QLfi) , 

q=u,d,s,c,b 

04 = (dLalfih/s) J2 i^LpYqia), 

q=u,d,s,c,b 

05 = {dlalpbLa) {qRpl^qRp), 

q=u,d,s,c,b 

q=u,d,s,c,b 

g _ 

0^ = T7r^daa^,y{mbR + mdL)baJ^'"', 

lOTT^ 

Os = -^^d^T^pa^,{7rn,R + m,L)hpG''^^ , 

lovr^ 

g _ _ 

= -^{dLalf,bLa){h^,l5l) , 
lOTT^ 

On = {dLal,,CLp){cRpl^hRa), 

O12 = {dLcl,iCLc){cR(3l^hRp), 

011 = {dLcl,j,UL(3){uRp'^^hRa), 

012 = {dLcl,iULa){uRpl^hRp) , (2) 

and the second operator set which are flipped chirality partners of the first: 

O'l = {dRa'^^CRp){cRp'^%Ra), 

02 = {dRa1iiCRa){cRi3YbRl3), 
O'l = {dRaJf,URf3){uRf3-f^bRa), 
O2 = idRalfj,URa)iuR(3l^bRf3), 

03 = (dRa-f^bRa) J2 (qR^i^qm^), 

q=u,d,s,c,b 

04 = (dRa'Jf.bRp) Y {qRHl^qRa), 

q=u,d,s,c,b 
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q=u,d,s,c,b 
q=u,d,s,c,b 

g _ 

O'g = J^{,dRa-i^,bRo,){ll^,l) , 

g _ _ 

O'lo = -^^{dRa-i^bRa){l-f^-f^l) , 

O'li = {dRa1f,CRp){cLI3l^bLa) , 

0'i2 = {dRali,CRa){cLpl^bLp) , 

O'li = {dRa1f,URp){uLI3l^bLa) , 

0'i2 = {.dRal^,URa){uL|3YbL|3) , (3) 

where a and (3 are S*?/ (3) colour indices and J^^'^ and Q^'^ are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 

The initial values for the first set of operators (eq.(|^)) ^] are 

Cl^3!..6,ll,12("^W/) = ' 

Cl'\mw) = 1, 

r^M, \ _ -^-^ In ^ I ox + /X 

^""^^ - 4(x-l)^^'''^+ 24(x-l)3 ' 

^'"("^^^ = -4(^3TjI^""+ 8(x-l)3 ' 

SM ' - 1 



CiT(^^) = -^(5(x)-C(x)) , 
sm fc'p^ 

Cl^..6,ll,12("^W^) = , 

rrit v^^ Vtb 

rrit v^^ Vtb 

+ ~~~~ i^N,tt + ^^,tcT^) (^^,66 + ^N,sbT^)G2{y) , 



^td ^tb 



mt v^^ Vtb 



and for the second set of operators (eq. (^] 



CZ^'uimw) = 0, 
CZs,n,i2imw) = 0, 

rrit i/^^ Vtb 

I __}_(tD Wb j_ FD Ws\fFU j_FU ^cbxp/ X 

^td ^tb 
T^tTTT'b Vtd ^td ^tb 

1 - V* - - - Vfs 
C'^{mw) = —2 {^§,bdT77 + (^N,bb + ^N,sbTr'>^^(y'> ' 
"^t ^td ^tb 

Cfoimw) = —A^NM^+^NM)(^N,bb + ^N,s>^^)Liiy), (5) 
rrit v^^ Vtb 

where x = / and y = mf / m'jj± . In eqs. (H) and (|^) we used the redefinition 



e"'" = \/^e"'"- (6) 

Here the Wilson coefficients Cf *^(mt^) and C^{mw) denote the SM and the additional charged 
Higgs contributions respectively. The functions -B(x), C(x), -D(x), Fi(2)(?/), G'i(2)(y), Hi{y) and 
Li{y) are given in appendix B. Note that in the calculations we neglect the contributions due 
to the neutral Higgs bosons since their interactions include negligible Yukawa couplings (see 
for details). 

Finally, the initial values of the Wilson coefficients in the model III (eqs. (^and (^) are 



ri2HDM 
'-^1,3, ...6,11, 12 


{mw) 


= 0, 








Q2HDM 


{mw) 


= 1, 








q2HDM 


{mw) 


= Cf'^l 


[mw) 




[mw 


ri2HDM 


{mw) 


= cl'H 


[mw) 


+ c^^ 


[mw 


ri2HDM 


{mw) 


= 


[mw) 




[mw 


^2HDM 
<-^10 


{mw) 


— Oio 1 


[mw) 




[mw 



Cl,2,3,...6,ll,12("^Vl/) - , 
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Cfo™(m^) = C'^,'\m^) + C[^{m^) . (7) 

These initial values help us calculate the coefficients (7^2H_da/ (ji2Hdm lower scale 

as in the SM (|^ references therein). The /z scale dependence of the coefficients in the LLog 
approximation can be found in the literature 0. The operators O5, Oq, Oh, O"^, 

O12 and ( Og, Og, O'^^^, O'j^^, 0'^2 O'^^^) give contribution to the leading order matrix 
element of 6 ^ 57 and the magnetic moment type coefficient Cj^^ {^) {C'y^^ {^)) is redefined in 
the NDR scheme as: 

rrib rUb 

+ ^^if ""''(Z^) + ^^^^^ • (8) 

uif) nil) 

Since (02"^) produce dilepton via virtual photon, their Wilson coefficient C2(yu) (C*2(/i)) 
and the coefficients Ci(/i), C^i^^), C6(/i) (C'2(/u), C'^ifJ'), C'g(/i)) induced by the operator 
mixing, give contributions to Cg^^ (fi) {Cg^^ (fi)). In a more complete analysis, one has to take 
into account the long-distance (LD) contributions, produced by real uu, dd and cc intermediate 
states, i.e. p, uj and i = 1,...,6 (Table |1]). These effects can be taken into account 
by introducing a Breit-Wigner form of the resonance propogator and it gives an additional 
contribution to Cg^ijj,) @, (Cg^^ifi)). Finally the effective coefficients Cg^ijj,) p|, and 
Cg^^{fi) are defined in the NDR scheme as: 

(3Ci(/x) + C^ifx) + 3C3(/i) + C4(/i) + 3C5(m) + C6(p)) 
+ A„<^ s) —K 2^ 
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h{l, s) {4Cs{fi) + 4C4(/x) + 3C5(/i) + C6(/i)) 

^/z(0, S) (C3(/i) + 3C4(p)) + ^ (3C3(/i) + C4(/i) + 3C,ifi) + CM) , (9) 

6 



and 



{3CM + CM + 3C',{fi) + C'M + 3C^(/x) + CM) 
+ Xu<h{z,s) - 2^ 2—— — 

- ^/.(l, s) {4CM + AC'M + 3CM + CM) 

- h{o, s) {CM + 3CM)) + 1 i^cM + CM + 3CM + cM) ■ m 



2 

where z = ^ and s = In the above expression, 77(5) represents the one gluon correction to 
the matrix element O9 with = The functions r/(s), uj{s), h{z, s) and h{0, s) are given 
in appendix C. In eqs. @ and ([10|) , the phenomenological parameter k = 2.3 is chosen to be 
able to reproduce the correct value of the branching ratio Br{B J/ipX Xll) = Br{B — > 
J/ipX) Br{J/i) ^ Xll) 13 . 

In the derivations of p and uo meson resonance effects, we used the dependence of the 
coupling fy- through the expression [R2 



/vs(g') = /vs(0)(^i + -p^(P{.,(0) + Pv,(g^)^ , (11) 

where the coupling fy^ is defined as < 0|g7^g|V^(g^)|0 >= /yj(5'^)e^, -fV,(0) and Py^iS^) are the 
subtraction constants (Table (^). The function Py.{q^) is 

^v.(a = T7^(-4-f r + 4(l + 2r)(i^)i/Mrctan(-^)V2A ^ (^2) 
loTT^r V 3 r 1 — r / 

where r = q'^/Am'^ and is the mass of the quark which produces the meson. This expression 

is valid in the region < < Arrig. For the q^ values, q^ > 4m^, we use the assumption ^2 

fvM')=h^^\) (Table®). 

Finally, neglecting the down quark mass, the matrix element for h de^e~ decay is obtained 



as: 



M = -^^V,,V;,|(C'9^^^(^)J7.(l-75)& + C^r^^(/i)rf7.(l + 75)&) 

+ [CMd^^.il-l5)b + C[Mdl^.il + l5)b) 6^756 (13) 
- 2('c7^^^(/i)^J^a^,g'^(l + 75)& + C^"^^(/i)^rf^a^.gMl-75)&) e^e] . 
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Y 


m / (GeV^ 


r(^ij ]+]-) (GeV) 




d.uy { 


o.zo iU 




(J . VJOVJ 




^(3) 


3.770 


2.6410"^ 


^(4) 


4.040 


7.2810-^ 


^(5) 


4.160 


7.80 10-^ 


^(6) 


4.420 


4.73 10-^ 



Table 1: Masses of mesons and decay widths r('0 — > Z+Z ) used in the calculations. 





fvmGeV) 


Mml,)iGeV) 


PviO) 


pm 


p 


0.162 


0.17 


-0.7498 


-0.0430 




0.166 


0.180 


-0.7744 


-0.0430 



Table 2: The decay couphngs and the substraction constants for p and u; mesons. 

3 The exclusive B 7re+e~ decay 
3.1 The formulation 

Now, we continue to present the differential decay rate and CP violating asymmetry in the 
process B — > 7re^e~ . To calculate the decay width, branching ratio, etc., for the exclusive B 
7re+e^ decay, we need the matrix elements (tt (i7^(l±75)6 B), and (tt dia ni,q" {1 ± '-)^)h B). 
Using the parametrization 



< vr(p^|rf7^(l ± 75)6|fi(pB) > 



{'^PB-q),U{q^) + q,f-{q^) , 

-{{2pB - q)^q^ - {ml - m2)g^}^;(g2) , (M) 



where pb and are four momentum vectors of B and tt mesons respectively and q = Pb — P-k, 
we get the double differential decay rate: 

dr{B vre+e-) Glal^m% \Vt,V;/ AV^^ 



d^/s dz 



2i%5 



Here 



2\i2 



(15) 



z') (16) 



and z — cosO , 6 is the angle between the momentum of the electron and that of B meson in 
the center of mass frame of the lepton pair, 



A = 1 + ^2 + ^2 _ 2^ _ 2s _ 2is , 



(17) 



2 2 

where t = — S- and s = 

For the form factors f+{q^) and v{q'^), we use the results due to the dispersion formulation 
of the light-cone constituent quark model p^] 



m-4 

1+ 



(1 _ 4)2.31 ^ ' 

where /+(0) = 0.24, i;(0) = 0.05 and m^^ = 6.71 GeV , m„ = 6.68 Gel/. 
Let us now turn to the CP-violating asymmetry, which is defined as 

dr(B->7re+e-) _ dV{B^Tie+e-) 

A _ _<iy/s dy^ /^qN 

'^^ dr{B~,TTe+e-) I dr(B^7fe+e-) ' 
d^/s dy/s 

The wilson coefficient C^^^ is the origin of the CP violating asymmetry since it is a function of 
Am = y^^y* ■ With the parametrization 



and using eq. (11^ we get 



C'i^^ = ^[ + \ui'2, (20) 



Acp = -2Im{K) ^ A (21) 

where 

A. = |/m(ere2)/+(g') + 2m,Im{Q{Cff + C^^^^} ^(g') !/+(?') I (22) 

and 

i\ = ei + ei 

^2 = 6 + ^2 (23) 

In our numerical analysis we used the input values given in Table (|^). 
3.2 Discussion 

In this section, we would like to study the dependencies of the differential Br, and Acp of 
the decay B ne^e^, for the selected parameters of the model III {^^tty ^mb) 5 using the 
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L CLl CLliiC uCi 


\I fi In 
V dl Lie 


■rric 




rrib 


/I fi ( r<\T\ 


-1 


1 on 


At 


U.U4 


i tot\Dd) 


o.yo ■ iU (^LreV ) 


rriBa 


D.Zo (^LreV J 


Trip 


U. ^ Oo (^LreV ) 




n 1 Qo ^'noA^^ 
U.ioy (^(orev j 




I to [KjQy j 




an 9fi iTIpVi 

OU.ZiU I VJC V I 




91.19 (GeV) 


Aqcd 


0.214 (GeV) 




0.117 


sinOy/ 


VO.2325 



Table 3: The values of the input parameters used in the numerical calculations. 



constraints |^ coming from the AF = 2{F = K,D, B) mixing ,the p parameter and the 
measurement by CLEO collaboration 



Br{B ^ X,7) = (2.32 ± 0.07 ± 0.35) 10" 



(24) 



In the calculations, we take ^^tc « ^mty ^mb ^Nij ~ where i or j are first or second 



generation indices (see |^ for details). Under this assumption the Wilson coefficients Cj, Cg 
and C[q can be neglected compared to unprimed ones and the neutral Higgs contributions are 
suppressed. 

In figs. |l| and ^ we plot the differential Br of the decay B ne^e^ with respect to the 
dilepton mass for the fixed values of ^§ = 40 and charged Higgs mass mjj± = 400 GeV 
at the scale fx = nii,. Fig. □ represents the case where the ratio \rth\ = \-^^\ « 1. It is shown 
that the differential Br obtained in the model III is smaller compared to the one calculated in 
the SM. Fig. ^ (|]) devoted to the case where r^^, >> 1 for the fixed value of f^^,;,, ^§bb — ^0 mi, 
i^Nbb — 90 mf,). The differential Br in the model III increases at this region (r^f, >> 1) and it 
enhances strongly compared to the SM with the increasing 



Now we present the values of Br for the B 
LD effects. After integrating over g^, we get 



Fig. I). 

TTe~^e~ decay in the SM and model III, without 



Br{B 



ne e 



0.62 X 10-*^ (SM) 



(25) 
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and for the model III 

f 0.27x10-^ (|r,,| «l,e^^,, = 40m,) 

0.54 X 10-7 {rtb » 1 , ^§,bb = 40 nib) (26) 

2.65 X 10-7 (rt, » 1 , = 90 m,) . 

Here, the strong enhancement of the Br can be observed for rtb » 1, especially with increasing 
^§iyf,- Note that, in the calculations of Br and the differential Br, we used the Wolfenstein 
parameters, p = —0.07, rj = 0.34. 

Figs. ^ and ^ show the dependence of Acp for the Wolfenstein parameters p = —0.07, rj = 
0.34, fixed values of ^^j,^ = 40 rrib and charged Higgs mass mH± = 400 GeV at the scale p = mb, 
for \rtb\ << 1 and >> 1 respectively. The CP violation in the model III for \rtb\ << 1 and 
C,Nbb — 4:0 nib is slightly greater than the one in the SM. However, it decreases for rtb » 1 
and becomes extremely smaller compared to the one calculated in the SM with increasing ^^bb 
(Fig. I). 

We also present < Acp > for two different Wolfenstein parameters in two different dilepton 
mass regions 



(p, v) 


SM 


model III 
^bb = 40mb 

\rtb\ « 1 


model III 
4^ = 40mb 
rtb » 1 


model III 
^bb = QO,mb 
rtb » 1 


q"^ regions 


(0.3,0.34) 


2.2010-^ 


2.2110-^ 


1.5810-^ 


0.7210-^ 


I 




0.6310-^ 


0.6310-^ 


0.4810-^ 


0.2410-^ 


II 


(-0.07,0.34) 


0.9910-^ 


1.1810-^ 


0.8210-^ 


0.3610"^ 


I 




0.3210-^ 


0.3210-^ 


0.2410"^ 


0.1110-^ 


II 



Table 4: The average asymmetry < Acp > for regions I ( 1 GeV < < mj/^ — 20MeV ) 
and II (mj/^ + 20 MeV < < m^/ - 20 MeV ) 

In conclusion, we analyse the dependencies of the differential Br, Br, Aqp and the average 
CP-asymmetry < Acp > on the selected model III parameters ( ^N,bby ^N,tt ) the decay 
B vre^e". We obtain that the strong enhancement of the differential Br (Br) is possible in 
the framework of the model HI and observe that Acp is sensitive to the model HI parameters 

i^N,bby ^N,tt)- 

4 The exclusive B pe+e~ decay 
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4.1 The formulation 



In this section ,we analyse the differential decay rate, Acp and ApB in the process B pe^e~ . 
At this stage, we need the matrix elements (p (i7^(l ±75)^ -B), and (p dia^^q^il ±75)6 B). 
Using the parametrization of the form factors as in the matrix element of the B — > pe+e~ 
decay is obtained as |l36| : 



M 



Gap„ 



+ 17^75^ [2Ctoie^^p,e*>^g" + iD^tot^l - iD^toti^* q){pB +Pp)p- iD^^toti.^* q)q^ 



(27) 



where e*^ is the polarization vector of p meson, pb and Pp are four momentum vectors of B 
and p mesons, q = Pb — Pp and 



Here 



Bi tot 

B2 tot 
-S3 tot 
Ctot 

Ditot 

D2tot 

Dstot 



A + A' , 
Bi + B[, 
B2 + B'2 , 

53 + ^3 , 

c + c" , 

Di + Z^'i , 
D2 + D'2, 
Ds + D' 



{2t 



A 

Bi 

B2 

Bs 
C 
Di 
D2 



q^ 

-C^,"f{q') + 2C,^^^^ {{ml - ml)g^{q') + q'g_{q')) 



q^ 



-Ci,g{q^) , 

-Cio/(g') , 
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Cio a-(g^) 



(29) 



and 



C 
D' 



-C[,9{q') , 
-C[oa+{q^) , 



(30) 



For the formfactors g{q^), a_{q^), a+{q^), g+{q^), g-iq"^), h^q"^), and /(g^) we use the dispersion 
formulation of the hght-cone constituent quark model |20| in the following pole form 



a^iq') 



0.036 



1 - 



q- 



2.75 ' 



(6.55)2^ 
0.03 



q^ 



2.85 ' 



(6.88)2^ 
0.18 



q^ 



2.73 ' 



-0.026 



3.04 ' 



(7.29)2^ 
-0.20 ' 



q^ 



2.76 ' 



(6.57)2^ 
0.003 



q^ 



3.42 ' 



(6.50)2 



(6.43)2 



1.10 



1 



q^ 



+ (- 



(31) 



Using eq. 



^ (5.59)2 ^(7.10)2^ ) 
we get the double differential decay rate: 



dq^dz 



212^5^^ 



2Xm% 



m%s{l + z')[\A 



Hot\ 



tot\ 



+ 



B 



2r 



Xml{l-z^) [\B2tot\^ + \D; 



2 tot 
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+ 



2r 



m 



B 



Itot 



2Xm%{l - r - - z') {Re {B.^^tBl^J + Re {D,t,tD*,J} 



{Re{B,totC:J + Re{AotDl,J} 



(32) 



where z = cos6 , 9 is the angle between the momentum of electron e and that of B meson in the 



center of mass frame of the lepton pair, X = 1 + t + s — 2t — 2s — 2ts, t = and s 



m 



B 



m. 



We continue to present the CP-violating asymmetry, which is defined as in eq. ([T9|) with 
the replacement of vr p. Using the same parametrization as in eq. (^) we get 

A. 



AcP = -2Im{Xu)-^X 



(33) 



where 



A. 



X , moA 



t 



+ 



.,.,,te^)|4.m|/(g^) + :^(6. + ^) + -^^^^ ' ^) ...2^ . 



s 



-/m(6) -4- 



X 



m,sg{q') g^{q') - ^ ((1 " t)g+{q') + s g-{q')) (1 + v^) 
1 - t - 



+ 



2mB 
rrib 
2mBt 



TUB X a+(g^) + ^^^(1 -t-s) 



TUB 



(34) 



and 



^ X{4mls{\Aot\' + \C\l,) + ^{6s + ^J{\Buot\' + \Duot\'' 



2r 



+ A^|(|i?2to*P + \D2t0t?) - X ^ ^ ' Re{BuotBl,,t + A*ot/^LJ} • (35) 

Finally, we present ApB which can give more precise information about the Wilson coefficients 
Cj^^ , Cg^^ and Cm. It is defined as: 



AMq') 



1 , dV 

dz 



, dV 
dz- 



dq^dz J -I dq^dz 
dz . „ . — h / (22;- 



dq^dz J-i dq^dz 



(36) 



After the standard calculation, we get 



^ 'n I 



(1 + Vi) g{q') ((1 - t)g^{q') + sg^{q')) + g^{q') f{q') 



l + t 



mB{l + Vi), 



(37) 
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4.2 Discussion 



In this section, we study the dependencies of the differential Br, Acp and ApB of the decay 
B pe^e^ for the selected parameters of the model III {^^tt^ ^Nbb)- the calculations, we 
use the same restrictions for the model III parameters, (see section 3) 

In figs. d^) we plot the differential Br of the decay B pe^e^ with respect to the 
dilepton mass for the fixed values of C,§f,h = 40 mb and charged Higgs mass mH± = 400 GeV 
at the scale p = nib, for the ratio \rtb\ = ItS^I << 1 {rtb = -rrr^ >> 1). The differential Br, 



^N,bb 



obtained in the model III, is smaller compared to the one calculated in the SM, for \rtb\ « 1. 
However, it increases at the region rtb » I and enhances strongly compared to the SM with 
the increasing ^^bb (Fig- 0); similar to the decay B vre+e^. To be complete, we present the 
values of Br for the B —>■ pe^e~ decay in the SM and model III, without the LD effects. After 



integrating over q , we get 



Br{B 



pe^e 



0.91 X 10~^ {SM) 



(St 



and for the model III 



Br{B — > pe~^e 



0.44 X 10-^ 
1.5 X 10-^ 



{\rtb\ « 1, ^_ 
{rtb » 1 , ^ 



D 

N,bb 



D 

NM 



3.2 X 10 



-7 



{rtb » 1 , 



NM 



= 40 rrib) 
40 nib) 
90 nib) ■ 



(39) 



The strong enhancement of the Br is observed for rtb » 1, especially with increasing ^^bb- 
Note that we used the Wolfenstein parameters, p = — 0.07,// = 0.34, in the calculation of Br 
and differential Br. 

Figs. |l^ and |TT| show the dependence of Acp for the Wolfenstein parameters p = 
—0.07, rj = 0.34, the fixed values of ^^bb — 4:0 rUb and charged Higgs mass mjj± = 400 GeV 
at the scale p = rUb, for \rtb\ << 1 and rtb » 1 respectively. The CP violation decreases in 
the region rtb » 1, especially with increasing ^§bb (Fig- 0)- Now, we give < Acp > for two 
different Wolfenstein parameters in two different dilepton mass regions 

Finally, we discuss ApB of the process under consideration. Figs. |l^ and [l^ show the 
dependence of App for the Wolfenstein parameters p = —0.07, rj = 0.34, the fixed values of 
^^bb = 40 m;, and charged Higgs mass = 400 GeV at the scale p = nib, for \rtb\ « 1 and 



rtb » 1 respectively. For rtb » 1 (Fig- llD ^fb changes its sign almost at s = 0.34, however 
for r^fe >> 1 (Fig. [1^) it is positive without LD effects. Therefore the determination of the sign 
of ApB in the region < s < 0.25 (here the upper limit corresponds to the value where ApB 
change sign in the SM) can give a unique information about the existence of the model HI. 
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C /\/f 


model III 
^bb — ^umb 
nb « 1 


model III 

^bb ~ ''''b 

\rtb\ > 1 


model III 

\rtb\ > 1 


q"^ regions 


(0.3,0.34) 


2.00 10"^ 


1.90 10"^ 


1.5010-^ 


0.5110"^ 


I 




0.60 10"^ 


0.57 10"^ 


0.5310"^ 


0.2510"^ 


II 


(-0.07,0.34) 


0.9710-^ 


1.00 10"^ 


0.7710-^ 


0.34 10"^ 


I 




0.3210"^ 


0.2910"^ 


0.2710"^ 


0.1410"^ 


II 



Table 5: The average asymmetry < Acp > for regions I ( 1 GeV < < mj/^ - 20 MeV ) 
and II (mj/^ + 20 MeV < < m^, - 20 MeV ). 

In conclusion, we analyse the selected model III parameters ( C^bb; ^N,tt ) dependencies 
of the differential Br ,Acp and Aps of the decay B pe+e~. We obtain that the strong 
enhancement of the differential Br is possible in the framework of the model III and observe 
that Acp and App are very sensitive to the model III parameters {^N,bb^ ^N,tt)- 

5 Conclusion 

We study the exclusive processes B 7ie^e~ and B pe^e~ which are induced by the 
inclusive h de^e~ decay. In such type of decays, it is informative to analyse the CP violating 
effects, in addition to the quantities like Br, ApB- The origin of the CP violation in the SM 
is the parameter = ^y^^* ■ In the model III, the couplings E,Yj and [| can also create the 
CP violation in case they are not real. However, in our work disregard this possibility not to 
enlarge the number of free parameters and we assume that the only CP violating effect comes 
from the CKM matrix elements, similar to the SM. 

Now, we would like to summarize the main results of our analysis: 

• The Br of the exclusive decays B 7re+e^ and B pe^e^ are sensitive to the model 
III parameters. In the region r^f, >> 1, a strong enhancement of the Br is observed with 
increasing in both decays eqs. As an example, BrnoLoiModel HI) ~ SBrnoLoiSM) 
for = 90 m;, ( P^lpSj ) and (^,^). Therefore their experimental investigations are a 
crucial test for the physics beyond the SM. 

• We calculated < Acp > for two different invariant mass region (see Table (§) and (^). 
We observe that < Acp > decreases with increasing for rtb » I in both regions. For 

"'^Here «, j and fc, I denote up and down quarks respectively. 
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rtb » 1 and — 90 m;,, < Acp > is rather smaller compared the one in the SM, for 
both regions (region I and II), i.e. < Acp >modeiiii ~ %30 < Acp >sm- 

• We calculated ApB for the B pe+e^ decay and observe that it does not change sign 
in the model III if the LD effects are excluded. This shows that the determination of the 
sign of ApB in the region < s < 0.25 will be informative to see the effects of the model 
III, if it exists. 

As a conclusion, the experimental investigation of the quantities we present here, will be an 
efficient tool to search for new physics beyond the SM. 
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Appendix 



A The essential points of the model III. 

The Yukawa interaction in the general case of 2HDM (Model III) is 



(40) 



where L and R denote chiral projections L{R) = 1/2(1 =p 75), (pi for i = 1,2, are the two scalar 
doublets, rjfj^ and ^fj^ are the matrices of the Yukawa couplings. The Flavor Changing (FC) 
part of the interaction can be written as 



^Y,FC = ^nQiL^UiR + ^nQiLhDjR + h.c. 



(41) 



with the choice of 0i and 02 



V2 







+ 



02 



_]_( V2H+ 
v/2 1, Hi + tH2 



(42) 



Here the vacuum expectation values are, 



<0,>=i=(O j ;<02>=O 



(43) 



and the couplings ^^'^ for the FC charged interactions are 



^ch ~ ^neutral VcKM , 
^ch — ^CKM ^neutral , 



(44) 



where ^n^utrai is defined by the expression 



(45) 



Note that the charged couplings appear as a linear combinations of neutral couplings multiplied 



by VcKM matrix elements (more details see 



^In all next discussion we denote £,^^^^^^11 ^^'^ ■ 
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B The necessary functions appear in the calculation of 
the Wilson coefficients 

The functions B{x), C{x), D{x), Fi(2)(?/), G'i(2)(y), Hi{y) and Liiij) are given as 



B{x) 
C{x) 
D{x) 

Fi{y) 
F2{y) 
Gi{y) 
G2{y) 
Hiiv) 



Li{y) 



—X 



+ 



X 



In X 



x — 1 (x — 1)2 
x/2-3 3x/2 + l 
X — 1 (x — 1)2 
19x736 + 25x736 -x76 + 5x73 - 3x'^ + 16x/9 - 4/9 



Inx 



+ 



(X- 1)3 

y{7-5y-8y'') y\3y-2) 
72(y-l)3 + 12(1/ -1)4 
y{5y-3) , y{-3y + 2) 



(x - 1)^ 



Inx , 



\ny 



+ 



12(y-l)2 Q{y-ir 
yi-y^ + 5y + 2) 



+ 



24(y - 1)3 4(y - 1)4 

y{y-3) , y 



\ny 



4(y-l)2 2(2/ -1)3 
1 — 4sm2^^(y xy 



In 2/ 



w 



y-l {y-lf 
A7y^ - 79y + 38 _ 3y^ - Qy + 4 
18(2/ - 1)^ 
1 



In 2/ 



108(2/ -1)= 



In 2/ 



i X2/ 
siv?6w 8 



+ 



2/-1 (2/-1)' 



In 2/ 



(46) 



C The functions which appear in the Wilson coefficients 

Cf and C^^^ 



The function which represents the one gluon correction to the matrix element Og is 

r/(s) = l + ^c.(l), 



TT 



(47) 



and 



2 2 4^ 2_, , 5 + 4S 

uj[s) = — TT i^Ms) msmll — s) ; 

^ ' 9 3 ^ ^ 3 ^ ^ 3(1 + 2s) 

2s(l + s)(l-2s) , . 5 + 9s-6s2 
m s + 



ln(l - s) - 



3(1 - s)2(l + 2s) 



6(1 - s)(l + 2s) 



(45 
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h{z, s) arises from the one loop contributions of the four quark operators Oi, Oq {0[, O'q) 



h{z, s) 



8, mb 8, 8 4 

— m mz -\ — - H — X 

9 /X 9 27 9 



2. ^, (la -in), iorx = ^<l 



2 arctan -t^^ , 



for X = ^ > 1, 



h{0,s) 



8 8, mfe 4, . 4. 

— m m s H — ITT , 

27 9 /X 9 9 ' 



(49) 



(50) 



where z — and s — 
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lO"? dBr/ds 



1.5 



1.25 





Figure 1: Differential Br as a function of q for fixed i^j,}, = 40 m;, in tlie region \rtb\ « 1, at 
the scale ii — nib for the process B — > 7re+e~. Here sohd hne and corresponds to the model III 
with LD effects, dashed line to the model III withouth LD effects and dotted dashed line to 
the SM withouth LD effects. 



10"7 dBr/ds 




Figure 2: The same as Fig 1, but at the region rtb » 1- 
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Figure 3: The same as Fig 2, but for fixed ^§^bb = 90 m;, value. 



. 1 



. 05 



-0 . 05 



-0.1 




Figure 4: Acp as a function of for fixed ^^j,^ = 40 nib in the region |rtb| << 1, at the scale 
/i = mb, for the process B — > ne'^e^. Here solid line corresponds to the model III with LD 
effects, dashed line to the SM withouth LD effects and dotted dashed line to the SM with LD 
effects. 
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Acp 




Figure 6: The same as Fig but for fixed ^^^^ = 90 nib value. 



10"7 dBr/ds 




Figure 7: Differential Br as a function of for fixed ^^hh = 40 mi, in the region \rtb\ << 1, at 
the scale n = rrib for the process B pe^e~ . Here solid line and corresponds to the model III 
with LD effects, dashed line to the model III withouth LD effects and dotted dashed line to 
the SM withouth LD effects. 



10*7 dBr/ds 




Figure 8: The same as Fig. ^ but at the region rtb » 1. 
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Figure 9: The same as Fig. ||, but for fixed ^^^f, = dOnif, value 



Acp 




Figure 10: Acp as a function of for fixed ^^1,1, = 4:0 nib in the region \rth\ « 1, at the scale 
fi = TTif,, for the process B pe^e" . Here solid line corresponds to the model III with LD 
effects, dashed line to the SM withouth LD effects and dotted dashed line to the SM with LD 
effects. 



Acp 





Figure 11: The same as Fig but at the region r^fe >> 1 
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Acp 




Figure 12: The same as Fig|Tl|, but for fixed (.j^hi, = 90 nib value 




Figure 13: ApB ^ function of for fixed ^^^b — 40 m;, in the region \rth\ « 1, at the scale 
/i = rrih for the process B — > pe^e~ . Here solid line and corresponds to the model III with LD 
effects, dashed line to the SM withouth LD effects and dotted dashed line to the SM with LD 
effects. 




Figure 14: The same as Fig. |13|, but at the region rtb » 1. 
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